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Abstract 

We performed density functional theory (DFT) 
calculations for a bi-layered heterostructure com¬ 
bining a graphene layer with a M0S2 layer with 
and without intercalated Li atoms. Our calcu¬ 
lations demonstrate the importance of the van 
der Waals (vdW) interaction, which is crucial for 
forming stable bonding between the layers. Our 
DFT calculation correctly reproduces the linear 
dispersion, or Dirac cone, feature at the Fermi 
energy for the isolated graphene monolayer and 
the band gap for the M0S2 monolayer. For 
the combined graphene/MoS2 bi-layer, we ob¬ 
serve interesting electronic structure and density 
of states (DOS) characteristics near the Fermi en¬ 
ergy, showing both the gap like features of the 
M0S2 layer and in-gap states with linear disper¬ 
sion contributed mostly by the graphene layer. 
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Our calculated total density of states (DOS) in 
this vdW heterostructure reveals that the graphene 
layer significantly contributes to pinning the Fermi 
energy at the center of the band gap of M0S2. 
We also find that intercalating Li ions in between 
the layers of the graphene/MoS2 heterostructure 
enhances the binding energy through orbital hy¬ 
bridizations between cations (Li adatoms) and an¬ 
ions (graphene and M0S2 monolayers). Moreover, 
we calculate the dielectric function of the Li inter¬ 
calated graphene/MoS2 heterostructure, the imag¬ 
inary component of which can be directly com¬ 
pared with experimental measurements of optical 
conductivity in order to validate our theoretical 
prediction. We observe sharp features in the imag¬ 
inary component of the dielectric function, which 
shows the presence of a Drude peak in the optical 
conductivity, and therefore metallicity in the lithi- 
ated graphene/MoS2 hetero structure. 

Introduction 

Currently, Li-ion batteries (LIB) are a very popular 
power source for portable electronic devices and 
hybrid vehicles. However, despite being renew¬ 
able and environmentally friendly, conventional 
Li-ion batteries are still far from meeting the na¬ 
tional US ABC target for a cheap ($ 150/KWh) 
and high capacity (3000 KWh/Kg) power source.- 
This is mostly due to the widespread use of 
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Figure 1 : Image of (a) graphene monolayer, (b) M0S2 
monolayer, (c) unit cell for graphene/MoS2 bilayer, and 
(d) unit cell for Li intercalated graphene/MoS2 bilayer; 
For a particular Li atom LU, the nearest Mo and C 
atoms are denoted with Moa and Ca; 


graphite based anodes in LIB, which only have a 
small theoretical specific capacity (372 mAh/g ).- 2 
Thus, the development of the next generation 
LIB calls for identifying alternative materials for 
the negative and positive electrodes. In this re¬ 
spect, both graphene and M0S2 monolayers have 
demonstrated potential as materials for design¬ 
ing superior electrodes with high specific capac¬ 
ity. High capacities from 600 to 1000 mAh/g 
has been observed for single atomically thin 
graphene nanosheets . 3 6 Moreover, several experi¬ 
mental works-- on transition metals and metal ox¬ 
ides supported on graphene as negative electrodes 
have shown enhanced performance for LIB. On 
the other hand, the usage of a transition metal sul¬ 
fide M0S2, with graphene like structure, was pro¬ 
posed as a potential anode material in a patent in 
1980 .- Despite its high insertion capacity 1000 
mAh/g), M0S2 alone has not been realized as the 
primary anode material due to its low cyclic stabil- 
ity.- 

In this Letter, we construct van der Waals 
layers— as the ‘intercalating compound ’ 4 using 
graphene and M0S2 bilayer. Our density func¬ 
tional theory (DFT) calculations show orbital hy¬ 
bridization between the layers, where the vdW 
interaction dominantly contributes to the stabil¬ 
ity of the heterostructure with significant binding 
energy. A shift in chemical potential is also ob¬ 


served in our calculations for the graphene/MoS 2 
bilayer when compared with the isolated graphene 
or M0S2 mono layers. Particularly noticeable is 
the linear dispersion in the DOS near the Fermi 
energy of the graphene/MoS2 bilayer, which re¬ 
sembles that of the isolated graphene monolayer. 
Thus, combining graphene with a M0S2 layer 
helps to pin the chemical potential near the gap 
center of M0S2 and creates semi-metal like elec¬ 
tronic structure in the heterostructure. 

Our DFT calculation shows an even larger bind¬ 
ing energy when the graphene/MoS2 bilayer are 
intercalated with Li ions. The higher stability in 
such heterostructure makes it a potential candidate 
for designing LIB anodes with larger specific ca¬ 
pacity and longer cycle life-time. The hybridiza¬ 
tion between carbon p, lithium s, and molyb¬ 
denum d orbitals causes the chemical potential 
to shift into the < 7 -bands making the Li interca¬ 
lated heterostructure fully metallic. The result¬ 
ing enhancement of the conductivity in lithiated 
graphene/MoS2 hetero structure is highly desirable 
for the anode materials in LIB. 

Earlier theoretical works— have demonstrated 
similar chemical potential shift using various 
adatoms in addition to Li, for single and bilayer 
graphene. For a related graphene and M0S2 het¬ 
erostructure, Miwa et al . 12 reported the DFT cal¬ 
culated density of states without observing any 
Fermi energy pinning or linear dispersion near 
the gap center of M0S2. In this Letter, we have 
demonstrated, for the first time, the interesting 
Dirac dispersion characteristics in the electronic 
DOS of graphene/MoS2 bilayer at the Fermi en¬ 
ergy. Although the lithiation makes the compound 
metallic as expected, discharging of LIB gradually 
turns the anode into semimetal while retaining the 
stability. In the language of electrochemistry, such 
processes can be presented by the following two 
half-reactions: - 3 

In anode: 


charge 

xLiCf, / M0S2 xLi + -\-xe~ -\-C(>/M0S2, 


discharge 

and in cathode: 

Li 1 -xCoOi + xLi 1 +xe 


charge 


LiCoOi 


discharge 


Therefore, our proof-of-principles calculations 
demonstrate the Li intercalated graphene/MoS2 bi¬ 
layer as potential candidates for designing high 
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performance anode materials for the next gener¬ 
ation LIBs. 

Computational Methods 

All calculations in this paper are performed using 
the plane-wave pseudo-potential code VASP—~— 
under the generalized gradient approximation 
of Perdew, Burke, and Ernzerhof (PBE). 17 For 
atomic core-levels, we have used projected aug¬ 
mented wave (PAW) potential s 1 — 1 ^ treating the 
2s2p of C, 2s of Li, 4p5s4d of Mo, and 3s3p of S 
as the explicit valence electrons. For all calcula¬ 
tions, the total energy during electronic relaxation 
is converged to 10 6 eV while the force/atom dur¬ 
ing ionic relaxation is converged to O.OleV/A. A 
maximum energy cutoff of 500 eV is used for 
plane-wave basis set. 

In the xy-plane, 6x6 unit cells of graphene and 
5x5 unit cells of M 0 S 2 are placed in a super-cell 
with a=b= 16 A, c=40 A, a — /3=90°, and y=120°. 
In our calculations, the Li atoms are intercalated 
between the graphene and M 0 S 2 layers at the cen¬ 
ter of the graphene honeycomb sites. For different 
metal ions, this location has been shown to have 
the largest binding energies by Chan et a/.— We 
have placed 12 Li atoms homogeneously in our 
super-cell. This sparse distribution of Li atoms 
approximates the interaction of isolated Li atoms 
with graphene and M 0 S 2 . The distance between 
Li atoms are large enough (ss 7 A) to minimize 
the orbital overlap between neighboring Li atoms. 

To incorporate the vdW interaction between 
the graphene and M 0 S 2 layers, we have used 
optB86b-vdW functional where the exchange 
functionals were optimized for the correlation 
part. 20 Therefore, the LDA correlation part present 
in the PBE functional is removed by using the pa¬ 
rameter AGGAC = 0.000 in the input file in order 
to avoid double-counting. 

Calculations for the isolated 6x6 graphene, 5x5 
M 0 S 2 , and 6x6 graphene + Li adatoms were per¬ 
formed using the same sized super-cell. For ionic 
relaxation of graphene + Li , graphene + M 0 S 2 , 
and graphene + Li + M 0 S 2 systems, we have used 
the r point to sample the Brillouin zone, while 
for all other calculations, e.g., DOS and dielectric 
functions, we used equally distributed 80 k-points 


in the irreducible Brillouin zone. 

Results and Discussion 

To estimate the adsorption energy, we define, 

AL — Ea+Eb — Ea+b, (1) 

where E^+b is the total energy of the A+B com¬ 
posite system, and E^/Eb is the total energy of the 
isolated A/B constituent system. 

In Table 1, we present our DFT calculated total 
energy for graphene, M 0 S 2 , the graphene/MoS 2 
heterostructure, and the lithiated graphene/MoS 2 
heterostructure. Using Eq. (1), we calculate the 
binding energy between graphene and M 0 S 2 bi¬ 
layer with and without Li ions in between. We 
have 25 Mo atoms in our supercell, and scaled 
the total adsorption energy to obtain AE per Mo 
ion. We find the AE per Mo ion is enhanced by 
about three times when Li intercalation is consid¬ 
ered. Without Li atoms, the graphene and M 0 S 2 
layers form a physisorbed heterostructure with « 
0.1873 eV adsorption energy per Mo ion, where 
the most important contribution comes from the 
surface van der Waals interaction. The adsorp¬ 
tion energy significantly increases in the presence 
of intercalated Li ions, and these higher adsorp¬ 
tion energies (0.6 eV per Mo ion or 1.26 eV per Li 
ion) suggest a more stabile structure. At the same 
time, the small change in atomic positions during 
the ionic relaxation clearly indicates the process 
is dominantly van der Waals driven physisorption. 
Thus, the Li ions are easily separable and lithiation 
process can be reversed. This is a highly desirable 
criterion for LIBs during ionization-deionization 
process which enhances the cycle lifetime of bat¬ 
teries. 

The electronic structure of the isolated graphene 
and M 0 S 2 monolayers is changed near the Fermi 
energy when they are assembled to form a van der 
Waals heterostructure. In [2} we present the DOS 
for the graphene/MoS 2 bilayer systems with and 
without lithiation and compare with the DOS of 
the constituent monolayers. A graphene mono- 
layer is a two dimensional Dirac material with 
linear dispersion near the Fermi energy. 2 A’- 22 This 
property manifests by the cone like feature in the 
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Table 1: Binding Energy of graphene and M 0 S 2 interface with and without Li intercalation using 
DFT+vdW 



Graphene 

MoS 2 

Graphene + 

Li 

Graphene/ 
MoS 2 

Graphene / 
MoS 2 + 

Li 

Total Energy (eV) 

-553.25348 

.414.14277 

-571.97493 

-972.07776 

-1001.27726 

Binding 
Energy (eV) 




0.1873 
(per Mo ion) 

0.6064 
(per Mo ion) 

1.2633 
(per Li ion) 
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Figure 2: (color online) (a) Total 
DOS (black) for graphene; (b) To¬ 
tal DOS (blue) for M 0 S 2 ; (c) Total 
DOS (red) for graphene/MoS 2 het¬ 
erostructure, orbital projected d- 
DOS of Mo (blue), and 72 -DOS 
(black) for C; and (d) Total DOS 
(red), d-DOS (blue) of Mo in 
M 0 S 2 , p-DOS (black) of C in 
graphene, and .v-DOS (green) of 
Li in the lithiated graphene/MoS 2 
heterostructure; the Fermi energy 
is at 0 eV, shown by the vertical 
dashed line. Total DOS in the unit 
cell (for (c) and (d)) are scaled 
down by a factor of 50 for visual 
clarity. Inset in (c) and (d) shows 
DOS near the Fermi level. 


DOS of graphene at the Fermi energy as shown in 
Ha). 

On the other hand, M 0 S 2 monolayer is a direct 
gap semiconductor with « 1.8 eV band gap . 23 In 
our DFT calculations, this band gap is estimated 
to be 1.4 eV as shown in [2tb). Due to the unfilled 
J-orbitals in Mo atoms, the semi-local DFT func¬ 
tionals are inadequate, and one can incorporate 
many-body corrections such as GW— or hybrid 
functionals HSE 25 which accounts for the missing 
0.4 eV band gap in the conventional DFT. In the 
present proof-of-principles calculations, our find¬ 
ings should remain unaffected from such band gap 
correction, while incorporating these corrections 
in large supercell with vdW interaction would 


greatly increase the computational cost. There¬ 
fore, we neglect such band-gap corrections in our 
calculations. 

We plot the total DOS of graphene/MoS 2 in H c ) 
with solid red curve. To provide visual clarity, to¬ 
tal DOS is scaled down by a factor of 50 while 
comparing with the partial rf-DOS (blue) of Mo 
in M 0 S 2 and p-DOS (black) of C in graphene. 
We find the presence of both the Dirac-cone-like 
feature from graphene and a gap-like feature from 
M 0 S 2 at the Fermi energy of the composite bilayer 
of graphene/MoS 2 - An enlarged version of these 
symmetric features in the DOS at Ep is presented 
in the inset of He). This finding can be inter¬ 
preted as the Fermi energy pinning at the gap cen- 
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Figure 3 : (color online) Calcu¬ 
lated real (blue) and imaginary 
(red) part of the dielectric func¬ 
tion for (a) graphene, (b) M0S2, 

(c) graphene/MoS2 bilayer, and 

(d) Li intercalated graphene/MoS2 
heterostructure; Due to the pro¬ 
portionality between optical con¬ 
ductivity and imaginary dielectric 
function, such results can be di¬ 
rectly compared with experimental 
measurements. See the main text 
for details. 


ter of M 0 S 2 with the Dirac states contributed from 
the graphene monolayer. Such pinning of E p can¬ 
not be identified in the isolated M 0 S 2 monolayer 
where the exact location of the chemical potential 
is arbitrary inside the gap. The effect of lithiation 
on the DOS is shown inOd). Using the similar 
scaling, we compare the total DOS with that of the 
partial d, p, and 5 orbitals of Mo, C, and Li atoms. 
The Mo and C atoms are selected from the M 0 S 2 
and graphene layers, which are closest to a partic¬ 
ular Li atom. These atoms are labeled in Qld). 
However, we find the DOS is no longer symmetric 
in the presence of Li adatoms. Due to a strong hy¬ 
bridization between Li s, Mo d and C p orbitals as 
well as the electron doping effect of Li, the chem¬ 
ical potential is shifted towards the unoccupied d 
bands, and thus, composite system becomes metal¬ 
lic. This implies that the loss of electrons from Li 
to the graphene and M 0 S 2 layers leads to a strong 
Coulomb interaction between Li ions and charged 
layers. It explains our observation of an increased 
binding energy with Lithiation. 

Our calculations for the real and imaginary part 
of dielectric function are presented in 0 Imagi¬ 
nary dielectric function has a direct relation to the 
optical conductivity via,— 

CO 

Re (Tap (a) = 4 ^Im £ a p(oo). (2) 

While the conducting nature of graphene can 


be understood by the sharp rise of the peak at 
the zero frequency as shown in [3ja) (solid red 
curve), the gap is evident for the isolated M 0 S 2 
monolayer as shown in [3lb). When the bilayer 
of graphene/MoS 2 is formed, the total conductiv¬ 
ity, as can be seen in [3£c), appears to be a lin¬ 
ear combination of contributions from each indi¬ 
vidual mono layers. On the other hand, lithia¬ 
tion causes metallicity in graphene/MoS 2 bilayer, 
where a large Drude peak appears near the zero 
frequency t[3}d)). 

Conclusions 

In this Letter, adsorption of graphene on a M 0 S 2 
monolayer is studied using DLT-based first- 
principles theory with and without intercalated 
Li atoms. Lormation of a stable heterostructure 
between the layers is observed, where the major 
contribution to binding comes from the vdW in¬ 
teraction. our calculations show that the addition 
of the graphene layer pins the Lermi energy near 
the gap center of M 0 S 2 . Precise identification of 
the Lermi energy opens up the potential appli¬ 
cations for LIB anode material through the ma¬ 
nipulation of electronic structure of various con¬ 
stituent element combination with graphene/MoS 2 
bilayer. Lithium insertion or intercalation between 
graphene/MoS 2 bilayer is found to have a bind- 
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ing energy within the physisorption range in our 
calculations, which demonstrates their prospect as 
superior anode materials. Therefore, our proof- 
of-principles calculations of adsorption energy, 
structural stability, electronic structure, and opti¬ 
cal properties in graphene/MoS 2 heterostructure 
provides a basis for future experimental and theo¬ 
retical study towards designing the next generation 
anodes for Lithium ion batteries. 
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